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^ ; ABSTRACT 
O ; 

^ ' Evidence for non-thermal activity in clusters of galaxies is well established 

^ ■ from radio observations of synchrotron emission by relativistic electrons. New 

• windows in the Extreme Ultraviolet and Hard X-ray ranges have provided for 
' more powerful tools for the investigation of this phenomenon. Detection of hard 

^ ■ X-rays in the 20 to 100 keV range have been reported from several clusters of 

^ ■ galaxies, notably from Coma and others. Based on these earlier observations 

• we identified the relatively high redshift cluster 1E0657-56 (also known as RX 
^ ■ J0658-5557) as a good candidate for hard X-ray observations. This cluster, also 
00 ' known as the bullet cluster, has many other interesting and unusual features, 

• most notably that it is undergoing a merger, clearly visible in the X-ray images. 
O ■ Here we present results from a successful RXTE observations of this cluster. We 
1^ • summarize past observations and their theoretical interpretation which guided 

us in the selection process. We describe the new observations and present the 

2 ■ constraints we can set on the flux and spectrum of the hard X-rays. Finally we 

^ ■ discuss the constraints one can set on the characteristics of accelerated electrons 
■ 

^ ■ which produce the hard X-rays and the radio radiation. 

^ I Subject headings: Galaxy Clusters: individual (1E0657-56, RX J0658-5557) - 

' Particle Acceleration: X-rays 



1. INTRODUCTION 

The intra-cluster media (ICM) of several clusters of galaxies, in addition to the well 
studied thermal bremsstrahlung (TB) emission dominating in the ~ 2 — 10 keV soft X-ray 
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(SXR) region, show growing evidence for non-thermal activity. This activity was first ob- 
served in a form of diffuse radio radiation (classified cither as relic or halo sources). The 
first cluster with diffuse emission detected in the radio band was Coma, and recent system- 
atic searches have identified more than 40 clusters with halo or relic sources. In the case 
of Coma, the radio spectrum may be represented by a broken power law (Rephaeli 1979), 
or a power law with a rapid steepening (Thierbach et al. 2003) or with an exponential cut- 
off (Schlickeiser et al. 1987). There is httle doubt that this radiation is due to synchrotron 
emission by a population of relativistic electrons with similar spectra, however, from radio 
observations alone one cannot determine the energy of the electrons or the strength of the 
magnetic field. Additional observations or assumptions are required. Minimum total (parti- 
cles plus field) energy or equipartition assumptions imply magnetic field strength of S ~ ^G, 
in rough agreement with the Faraday rotation measurements {e.g.YJivn et al. 1990), and a 
population of relativistic electrons with Lorentz factor 7 ~ 10*^. In the papers cited above, 
it was also realized that because the energy density of the Cosmic Microwave Background 
(CMB) radiation (temperature Tq) mcmb = 4 x 10^^'^(To/2.8 K)'^ergcm~^ is larger than the 
magnetic energy density mb = 3 x 10~^''(i?//xG)^ erg cm~^, most of the energy of the rela- 
tivistic electrons will be radiated via inverse Compton (IC) scattering of the CMB photons, 
producing a broad photon spectrum (similar to that observed in the radio band) around 50 
keV (for 7 ~ 10^) . Thus, one expects a higher flux of non-thermal X-ray radiation than radio 
radiation. Detection of HXRs would then break the degeneracy and allow determination of 
the magnetic field and the energy of the radiating electrons. Moreover, since the redshift 
dependence of the CMB photons is known, then in principle, the cosmological evolution of 
these quantities can also be investigated. 

While the dctcctability of such radiation would be easiest in the soft X-ray range, where 
very sensitive imaging instruments are available, in reality, in this range, it is masked by the 
prominent thermal brcmsstrahlung emission, a general characteristic of clusters. With this, 
the most promising is either the hard X-ray (HXR) band, beyond the energies where TB fiux 
dominates, or alternatively, in the very soft X-ray or extreme ultraviolet regime. Recently 
HXR emission (in the 20 to 80 keV range) at significant levels above that expected from 
the thermal gas was detected by instruments on board Beppo^KX. and RXTE satellites 
from Coma (Rephaeli et al. 1999. Fusco-Femiano et al. 1999, Rephaeli & Gruber 2002, Fusco- 
Femiano et al.2004i), Abell 2319 (Gruber & Rephaeh 2002), Abell 2256 (Fusco-Femiano et 
al. 2000, Rephaeli & Gruber 2003, and Fusco-Femiano, Landi, & Orlandini 2005), and a 
marginal (~ 3cr) detection from Abell 754 and an upper limit on Abell 119 (Fusco-Femiano 



^The results of this paper have been chaUenged by an analysis performed with different software by 
Rossetti & Molendi (2004). 
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et al. 2003). We also note that a possible recent detection of non-thermal X-rays, albeit at 
lower energies, has been reported from a poor cluster IC 1262 by Hudson et al. (2003). 
All those clusters are in the redshift range 0.023 < z < 0.056. Notable recent exception 
at a higher redshift is Abell 2163 (Rephaeli, Gruber, & Arieli 2006) where the reported 
nonthermal flux is consistent with the upper limit set by BeppoSAX (Feretti et al. 2001). 
It should also be noted that excess radiation was detected in the 0.1 to 0.4 keV band by 
Rosat, BeppoSAX and XMM-Newton and in the EUV region (0.07 to 0.14 keV) similar excess 
radiation was detected by the Extreme Ultraviolet Explorer from Coma (Lieu et al. 1996) 
and some other clusters. A cooler (kT ~ 2 keV) component and IC scattering of CMB 
photons by lower energy (7 ~ 10^) electrons are two possible ways of producing this excess 
radiation. However, some of the observations and the emission process are still controversial 
(see Bowyer 2003). 

Here we present results from ~ 309 ks RXTE observations of cluster 1E0657-56 (also 
known as RX J0658-5557) with a considerably higher redshift of ^ = 0.296, which manifests 
many other interesting features (Markcvitch 2005). All those hard X-ray observations - 
especially with the addition of 1E0657-56 - encompass a wide range of temperature, redshift 
and luminosity, indicating that HXR emission may be a common property of all clusters 
with significant diffuse radio emission. In the next section we give brief descriptions of the 
emission processes and the considerations which led to the selection of 1E0657-56 as a target 
for HXR observations by RXTE . In §3 we describe the observations and the results from 
our spectral fits. Finally in §4 we discuss the significance of our results and their implication 
for the radiation and acceleration mechanisms. 



2. EMISSION PROCESSES AND TARGET SELECTION 

As stated above electrons of similar energies can be responsible for both the IC— HXR 
and synchrotron— radio emission and the ratio of these fiuxes depends primarily on the ratio 
of the photon (CMB in our case) to magnetic field energy densities. For a population of 
relativistic electrons with the spectrum A^(7) = ^7"^ = A^totai(p — l)7~^7mTn ; for7 > 7inin, 
the spectrum of the emitted luminosities of both radiation components is given by 

uLi{u) = CroA^total7m7n^iiA(p)(i^/i^cr,i)^~^, (1) 
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where ro is the classical electron radius, c is the speed of light and F = (p + l)/2 is the 
photon number spectral index^. 

For synchrotron i^cr.synch = '^I'Lin^B^/'^, with = eB±/{2Trmc) and ligynch = -BV(87r), 
and for IC, Ui is the soft photon energy density and i^cr,ic — 7min(^^)- For black body 
photons uic — (87r^/15)(A;T)^/(/ic)^ and (/ii/) = 2.S kT. Ai are some simple functions of the 
electron index p and are of the order of unity (see e.g'. Rybicki & Lightman 1979). Because 
we know the temperature of the CMB photons, from the observed ratio of fluxes we can 
determine the strength of the magnetic field. For Coma, this requires the volume averaged 
magnetic field to be i? ~ 0.1 /xG, while equipartition gives B ~ 0.4 /xG while Faraday rotation 
measurements give the average line-of-sight field of i?i ~ 3 /xG (Giovannini et al. 1993, Kim 
et al. 1990, Clarke et al. 2001; 2003). (In general the Faraday rotation measurements of 
most clusters give B > //G; see e.g. Govoni et al. 2003.) However, there are several factors 
which may resolve this discrepancy. Firstly, the last value assumes a chaotic magnetic field 
with scale of few kpc which is not a directly mcasTired quantity (see e.g. Carilli & Taylor 
2002). Secondly, the accuracy of these results have been questioned by Rudnick & BlundcU 
(2003) and defended by Govoni & Feretti (2004) and others. Thirdly, as pointed by Brunctti 
et al. (2001), a strong gradient in the magnetic field can reconcile the difference between 
the volume and line-of-sight averaged measurements. Finally, as pointed out by Petrosian 
2001 (POl for short), this discrepancy can be alleviated by a more reahstic electron spectral 
distribution (e.^-. the spectrum with exponential cutoff suggested by Schhckeiser et al. 1987) 
and/or a non- isotropic pitch angle distribution. In addition, for a population of clusters 
observational selection effects come into play and may favor Faraday rotation detection in 
high B clusters which will have a weaker IC flux relative to synchrotron. 

A second possibility is that the HXR radiation is produced via bremsstrahlung by a 
population of electrons with energies comparable and larger than the HXR photons. If a 
thermal distribution of electrons is the source of this radiation, such electrons must have a 
much higher temperature than the gas responsible for the SXR emission. For production 
of HXR fiux up to 50 keV this requires a gas with kT > 30 keV and (for Coma) with an 
emission measure about 10% of that of the SXR producing plasma. Heating and maintaining 
of the plasma to such high temperatures in view of rapid equilibration expected by classical 
Spitzer conduction is problematical. It has also been suggested by various authors (see, 
e.5'. EnBlin, Lieu, & Biermann 1999; Blasi 2000) that the HXR radiation is due to non- 
thermal bremsstrahlung (NTB) by a power law distribution of electrons in the 10 to 100 keV 



^ These expressions are valid for p > 3 or F > 2. For smaller indices an upper energy limit 7max also must 
be specified and the above expressions must be modified by other factors which are omitted here for the sake 
of simplification. 
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range. However, it was shown in POl that the NTB process faces a serious difficulty, which 
is hard to circumvent, because compared to Coulomb losses the bremsstrahlung yield is very 
small; Ibrem ~ 3 X 10~^(£'/25 keV)^/^ (see Petrosian 1973). Thus, for continuing production 
of a HXR luminosity of 4 x 10^^ erg s^^ (observed for Coma), a power of i^HXR/^brcm ~ 10^^ 
erg s~^ must be continuously fed into the ICM, increasing its temperature to T ~ 10^ K 
after 3 x 10'' yr, or to 10^° K in a Hubble time^. Therefore, the NTB emission phase must 
be very short lived. A possible way to circumvent the rapid cooling of the hotter plasma by 
conduction or rapid energy loss of the non-thermal particles is to physically separate these 
from the cooler ICM gas. Exactly how this can be done is difficult to determine but strong 
magnetic fields or turbulence may be able to produce such a situation. 

In what follows we shall assume the IC process as the working hypothesis for the origin 
of the HXR fiux. However, in view of the above mentioned difficulties more observations 
are acutely needed to determine the HXR emission process. Such observations are chal- 
lenging: thermal emission from clusters of galaxies can be well characterized as thermal 
bremsstrahlung with emission lines due to atomic transitions, with a characteristic temper- 
ature of a ~ a few to ~ 15 keV. This means that any reliable detection of the non-thermal 
component relies on an instrument sensitive above at least 10 keV. Sensitive measurements 
require instruments employing focusing optics, but no such instruments are currently oper- 
ational. We must thus rely on the coUimated detectors on-board of 5e^|)oSAX, RXTE, or 
Suzaku HXD, but for such instruments, good limiting sensitivity requires long observations. 
Therefore, the target selection - based on the good prediction of the HXR fiux - requires 
careful considerations of all aspects of the phenomenon. 

The most important selection criterion is the presence of diffuse radio emission. About 
40 clusters are known to have such emission, which can be classified as halo or relic (see 
Giovannini, Tordi, & Feretti 1999, Giovannini & Feretti 2000, Feretti et al.2000). The 
fraction of clusters with diffuse radio emission rises with the SXR luminosity (Giovannini 
et al. 1999) which is in turn correlated with temperature, via Lsxr oc (see e.g. Allen 
& Fabian 1998), indicating a correlation of the non-thermal activity with temperature. In 
addition, turbulence and shocks, present in merging clusters, are the most likely agents of 
acceleration and could be the cause of the higher temperatures, as well as the nonthermal 
activity. Consequently, high temperatures and presence of substructure also must infiuence 
the selection (see Buote 2001 and Schuecker et al. 2001). Note that higher temperatures 
make the detection of the HXR fiux more difficult. However, at higher redshifts this effect 



^ These estimates are based on energy losses of electrons in a cold plasma which is an excellent approxima- 
tion for electron energies E ^ kT . As E nears kT the rate of loss of energy decreases and the bremsstrahlung 
yield increases. For E/kT > 4 this increase will be at most about a factor of 2. 
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is offset by tlie spectral redsfiift z = Z — \. Based on tfie criteria given above, a list of some 
of the most promising clusters are given in Table 1. 



Table 1 

OBSERVED AND ESTIMATED PROPERTIES OF CLUSTERS 



Cluster 


z 


kT" 
keV 


-^1.4GHz 

mJy 


0C,b 

arcmin 


FsXR 

Fi 


HG 


-^HXR 


Coma 


0.023 


7.9 


52 


30 


33 


0.40 


1.4(1.6) 


A 2256 


0.058 


7.5 


400 


12 


5.1 


1.1 


1.8(1.0) 


1E0657-56 


0.296 


15.6 


78 


5 


3.9 


1.2 


0.52(0.5) 


A 2219 


0.226 


12.4 


81 


8 


2.4 


0.86 


1.0 


MACSJ0717 


0.550 


13 


220 


3 


3.5 


2.6 


0.76 


A 2163 


0.208 


13.8 


55 


6 


3.3 


0.97 


0.51 


A 2744 


0.308 


11.0 


38 


5 


0.76 


1.0 


0.41 


A 1914 


0.171 


10.7 


50 


4 


1.8 


1.3 


0.22 



" Prom Allen & Fabian (1998), except 1E0657-56 data from Liang et al. (2000) 
From Giovannini et al. (1999, 2000), except 1E0657-56 data from Liang et al. (2000) 

Approximate largest angular extent. 

Estimates based on equipartition. 
^ Estimates assuming C7min = 10®, with observed values in parentheses for Coma from 
Rephaeli et al. (1999; 2002), and Fusco-Femiano et al. (1999; 2004) and for Abell 2256 by 
Fusco-Femiano et al. (2000; 2005) and Rephaeh & Gruber (2003). 
^ Fq = 10~^^ erg cm"^ s~^ 

These are only qualitative criteria but for the IC model one can give some quantitative 
estimates. At a redshift z = Z — 1 the CMB energy -u oc and the critical frequency 
i^cic o: Z. As a result the ratio of IC to synchrotron fluxes 

n = FHXR/Fradio OC Z^+^ / B^, (2) 

SO that for an observed radio flux one would like to choose clusters with the lowest magnetic 
field B and the highest redshifts. Unfortunately the B field strengths in most of the diffuse 
radio emitting clusters are not known (Coma is an exception) so that one must rely on some 
theoretical arguments to estimate the value oi B. If we assume some proportional relation 
{e.g. equipartition) between the energies of the magnetic field and non-thermal electrons 
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C^e = A/total ^TminmC = Q- — , (3) 

p — 2 Stt 3 

where R = OdA/'i is the radius of the (assumed) spherical cluster with the measured an- 
gular diameter 6 and angular diameter distance dA{Z), and equipartition with electrons is 
equivalent to C = 1- 

From the three equations (1), (2) and (3) we can determine the three unknowns B,Se 
(or A^totai) and Fhxr purely in terms of C)7min, and the observed quantities (given in Table 
I) z,9 and the radio flux 



i^F^^dM = crlN,,,^iA,,^,^{p){p - l)7j;T;(5V87r)(Zz//z/er,.,)'-7(47r4(^)) (4) 

where = dAZ^ — Zr{Z) is the luminosity distance to a source with the co-moving 
coordinate r{Z). The result is^ 

65^-3/2^/lo3^^ feYf b Y fr{z)^^ 



and 

where we have defined Fq = 10 erg cm ^ s ^. Note also that in all these expressions one 
may use the the substitution R = 3.39 Mpc {O/b') r{Z)/Z. 

To obtain numerical estimates of the above quantities in addition to the observables 
-^radio) ^ and redshift z we need the values of ^ and 7min- Very little is known about these two 
parameters and how they may vary from cluster to cluster. From the radio observations at 
the lowest frequency we can set an upper limit on 7min; for Coma e.g. ^min < 4 x 10'^(/iG/i?). 
We also know that the cut off energy cannot be too low because that will require excessive 
amount of energy (for p > 3) which will go into heating of the ICM gas via Coulomb collisions. 



"^Here we have set the Hubble constant Hq = 70 km Mpc~^ s~^, the CMB temperature Tq = 2.8 K, and 
the radio frequency u = 1.4 GHz. In general B'^+^ oc Hoiy^~^ , N^otsd and Fhxr oc HqTq^^. We have 

also assumed an isotropic distribution of the electron pitch angles and set B = Sj_(4/7r). 
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A conservative estimate will be 7inin ~ 10^. Even less is known about (. The estimated values 
of the magnetic fields B for the simple case of F = 2, equipartition [i.e. C = 1) ^-nd low energy 
cut off 7niin = 10^ are given in the 7th column of Table 1. As expected these are of the order 
of a few fiG; for significantly stronger field, the predicted HXR fiuxcs will be below what 
is detected (or even potentially detectable). For F = 2 the magnetic field B oc (C7min)~^^^ 
and Fhxr oc (C7mm)^^^ so that for sub-/xG fields and Fhxr ~ -^o we need C7min ^ 10^. 
Assuming 7min = 10^ and C = 10^ we have calculated the expected fluxes integrated in the 
range of 20 - 100 keV (which for F = 2 is equal to 1.62 x (20keVFHXR(20keV)) shown on 
the last column of Table 1. The variation of this flux with redshift based on the observed 
parameters, 9 and -Fradio(i^ — 1-4 GHz) of Coma and A2256 are plotted in Figure 1 for three 
values of F = 1.75, 2.0 and 2.25 {p = 2.5, 3 and 3.5) and assuming a constant metric radius 
R which a reasonable assumption. We also plot the same assuming a constant angular 
diameter. This could be the case due to observational selection bias if most diffuse emission 
is near the resolution of the telescopes. These are clearly uncertain procedures and can give 
only semi-quantitative measures. However, the fact that the Coma and A2256 have the 
highest fluxes, and that our observations of 1E0657-56 described below, yield a flux close to 
our predicted value is encouraging. Clearly, all the other clusters in Table 1 are similarly 
promising candidates for future HXR observations. 



3. OBSERVATIONS AND ANALYSIS 
3.1. Description of Observations 

Cluster 1E0657-56 was observed with RXTE during 95 separate pointings from 2002 

August 19 to 2003 April 3, collecting approximately 400,000 seconds of data. After filtering 
the PCA data by following the standard selection criteria, we obtained a total of 309,630 
seconds of good PCA data. During the campaign, only PCU and PCU 2 were on. Since 
PCU had lost its propane guard layer, in order to ensure that we have the best possible 
calibration, we used only the PCU 2 detector. The selection criteria involved excluding 
data obtained when the detector pointing was more than 0.02° off the source, excluding data 
obtained when the Earth elevation angle was less than 10°, and excluding times near passage 
through the South Atlantic Anomaly which caused large variations in the count rate. 

The amount of good time for each of the two HEXTE detector clusters is about one- 
third that for the PCA data: 105,450 seconds for Cluster A and 105,510 seconds for Cluster 
B. The ratio of HEXTE net time to PCA net time is reasonable, since each HEXTE cluster 
spent half of the time for background measurement and a fraction of the time was lost due 
to electronic dead time caused by cosmic rays. 
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Dotted (black): Coma, Constant Ang. Dia. 
Long dashed (red): A2256, Constant Ang. Dia. 
Short dashed (red): A2256, Constant Radius 
Solid (black): COMA, Constant Radius /''^ 
=103, ^=103 >^ 

' min ' J 




0,01 



0.1 



Redshift 



Fig. 1. — Predicted variations of the HXR flux with redshift assuming a constant metric (solid and dashed) 
or angular (dotted and long-dashed) diameters, using COMA (black and) and A2256 (red) parameters 
assuming C = 7niin = 10'^. In each group the IC photon spectral index F — 2.25,2.0,1.75, from top to 
bottom. Filled circles based on observations and stars based on the estimates given in Table 1. 
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The PCA background was estimated using the ftool pcabackest (v 3.0) with the 
"L7/240" model provided by the RXTE GOF (Guest Observer Facihty). The estimator 
program used the model known as pca_bkgdcmf aintl7_eMv20051128.mdl as the background 
model, which in turn took into consideration the instantaneous particle-induced intrinsic 
instrument background and activation - mainly due the SAA (South Atlantic Anomaly) 
passages - as well as the the cosmic X-ray background. We note here that the above back- 
ground model, available from HEASARC at NASA's GSFC (Dr. Craig Markwardt, priv. 
comm.) allows robust estimate of background up to ~ 30 keV even for faint sources such as 
1ES0657-56. 

We used the top layer of PCU 2 in the 3 — 30 keV range. The mean background 
subtracted PCU2 count rate was 1.418 ±0.007 count s~^ in the top layer over this band. For 
the response matrix, we used the standard PCA response matrix, generated via the pcarsp 
software tool (ver. 10.1), as appropriate to the beginning of the observations (2002 August); 
while the RXTE PCA response matrix changes with time, the change is gradual, and the 
values of spectral fits obtained using the matrix appropriate for the end of the observations 
(2003 April) showed no discernible difference. For HEXTE, the background subtracted count 
rates over the 20 ~ 70 keV energy band were 0.0103 ± 0.0227 count s^^ and 0.0700 ± 0.0192 
count s~^ for HEXTE detector clusters A and B, respectively; above 70 keV, the source was 
clearly not detected. We used the standard redistribution files, xh97mar20c_pwa_64b . rmf 
and xh97inar20c_pwb013_64b.rmf , with the standard HEXTE effective area files available 
through HEASARC. 

While the RXTE data clearly detected the cluster, an inclusion of higher spectral 
resolution data allows a better constraint on the temperature of the thermal component and 
thus better constraints on the nature of the non-thermal emission. To this end, we considered 
the published ASCA data, but because the best statistical accuracy was provided in the 
XMM-Newton data, we extracted and analyzed archival spectra collected during the XMM- 
Newton pointing in 2000 October 20-21. We reduced the data using a procedure described 
in Andersson & Madejski (2004), but here, we used the XMMSAS_20050815 release of the 
XMM-Newton data analysis software. We screened the data against any obvious flares; this 
meant any segments of data with the total count rate greater than 12 counts s~^ for pn, and 
4.5 and counts s~^ for either MOS. After cleaning, the effective exposure times are 24,916 
s for the MOSl data, 24813 s for the M0S2 data, and 21,704 s for the pn data. We note 
that the detailed analysis of the XMM-Newton data for this cluster - and in particular, the 
analysis of the spatial structure - will be presented in Andersson et al (in preparation) . For 
all three XMM-Newton detectors, we extracted source counts from a region of 4' radius 
around the centroid of the count distribution. For background, we used an annulus with 
inner and outer radii of 8' and 12', with all obvious point sources removed. We also removed 
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the data corresponding to the Cu K hne in the pn camera, corresponding to 7.9 to 8.2 keV 

spectral range. The data were grouped to include at least 40 counts per bin for the pn 
data files, and 25 counts per bin for the MOS data files. We used the standard instrument 
resolution matrices and effective area files as released in the SAS version as above. In the 
subsequent analysis, we restricted the two MOS and pn data to the I - 10 keV range, to 
avoid any residual cross-calibration problems between the three instruments. 

3.2. Analysis and Spectral Fitting 

3.2.1. RXTE Data 

We first present an analysis of the RXTE observations and then a joint analysis of 
the RXTE and XMM-Newton data. In all subsequent fits we fix the redshift at 0.296. 
The models consist of a mekal thermal emission component, where we also allow a second 
component, either thermal, or a power law. We first use a single thermal mekal model to 
fit the RXTE data. The model includes a cold absorption due to neutral gas in the line 
of sight in our Galaxy. Since the RXTE data alone do not constraint the Galactic column 
density, we fix this parameter at 4.6 x 10^° cm~^, the value from radio measurements, but 
also consistent with the best fit to ASCA data (Liang et al. 2000). We note that all errors 
quoted correspond to 90% confidence limits. 

The best fit temperature under an assumption of a simple, one-temperature mekal model 
is 12.1 ±0.4 keV. This is in moderate agreement with the value 14.5l^ 7 keV determined from 
a joint fit to the ROSAT and ASCA data (Liang et al. 2000), and the value keV 
from Chandra (Markevitch et al. 2002). The metal abundance is A — 0.16 ± 0.04 Solar, 
which is moderately lower than the value A = 0.33 ± 0.16 from the joint ROSAT and ASCA 
data but is in good agreement with the Chandra value 0.11 ±0.11. The best-fit model yields 
the of 114 for 98 degrees of freedom. The ratio of the RXTE data to the best-fit mekal 
model indicates systematic upward rise above 10 keV, suggesting that an extra component 
may be required. If we adopt this secondary component to be a power-law, then is 102 
(96 degrees of freedom), with poorly determined index (F > 2). The 90% confidence interval 
of the power-law fiux in the 20-100 keV range is (0.3 ± 0.2) Fq. We thus conclude that even 
from the RXTE data alone, we have a marginal but suggestive detection of hard X-ray fiux 
in 1E06567-56. 



3.2.2. XMM-Newton Data 



The XMM-Newton data extracted as above were first fitted to a model including ab- 
sorption due to neutral gas with Solar abundances and a thermal (mekal) plasma. Since 
the absolute calibration of the three XMM-Newton detectors might vary, we allowed the 
respective normalizations not to be tied to each other. The best fit absorption was close to 
the value inferred from the radio data of Liang et al. (2000), with the column 2.8 ± 1.0 x 10^° 
cm~^, but here, inclusion of the MOS data below 1 keV altered the resulting best fit column 
to 5 ±0.5 X 10^° cm~^, We concluded that the column derived from the XMM-Newton data 
is consistent with 4.6 x 10^'^ cm~^ cited by Liang et al. (2000), and used that value in all 
subsequent fits. The temperature of the mekal component was 12.0 ±0.5 keV, and elemental 
abundances 0.24 ± 0.04 Solar, consistent with all data sets above. The fit is acceptable, at 
= 1367 for 1408 d.o.f. 



Table 2 

PARAMETERS FROM SPECTRAL FITTINGS 



Data Set 


Parameter 


Single Thermal 


Double-Thermal 


Thermal+Power Law 




(nn/lO^" cm-2) 


4.6^ 




4.6-^ 




kTi (keV) 


12.1 ±0.4 




11.7 ±0.5 




kT2 (keV) 








RXTE 


Abundance (Solar) 
Photon Index 

j^lOOfceV / TP 


0.16± 0.04 




0.25 ± 0.08 

> 2 
0.3 ±0.2 




xVdof 


114/98 




102/96 




(n^/lO^o cm-2) 


2.8 ± 1.0 


4.6^' 


4.6-'' 




kTi (keV) 


12.1 ±0.2 


10.1 ±0.9 


11.2 ±0.8 


RXTE and 


kT2 (keV) 




50 (> 30) 




XMM 


Abundance (Solar) 
Photon Index 


0.19 ± 0.03 


0.19 ± 0.03 


0.22 ± 0.04 
1.6 ±0.2 




rplOOkeV / rp 




0.5 ±0.3 


0.5 ±0.3 




xVdof 


1483/1508 


1471/1506 


1464/1506 



Fq = 10 erg cm ^ s ^ 

f denotes parameter fixed at the given value 
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3.2.3. Joint RXTE and XMM-Newton Data 

In the analysis below, we include both the RXTE and XMM-Newton data, since 
the very good effective area of the instrument coupled with good spectral resolution of its 
detectors allows tighter constraints on spectral parameters for the emission detected below 
~ 10 keV, and thus mainly on the thermal component. However, the only reliable approach 
here is to perform the spectral fitting simultaneously. 

The joint analysis of the RXTE and XMM-Newton data provide more evidence of 
the need for a secondary component. We use the XMM-Newton and RXTE data over the 
bandpasses as above; we allow the normalization of the RXTE instruments to be different 
than that for the XMM-Newton instruments, which in turn are allowed to vary among 
themselves. A single isothermal (mekal) fit yields an adequate fit, with of 1483 for 1508 
degrees of freedom. The hydrogen column density uh is (2.8 ± 1.0) x lO^'' cm~^, which is 
marginally lower than the Galactic value of Liang et al. (2000) (see the discussion above); in 
the subsequent fits we adopt the value of Liang et al. (2000) of 4.6 x 10^*^ cm^^. Regardless 
of the exact value of absorption, the plasma temperature is 12.1 ± 0.2 keV, and elemental 
abundances are 0.19 ± 0.03 Solar. We note that the true errors on those quantities are only 
approximate, since as pointed out by Markevitch et al. (2002; 2004), this cluster shows some 
temperature structure, while we use an average temperature. 

Adding a power-law component to the isothermal model improves the value, which is 
now 1464/1506 d.o.f. Such change in - with addition of two parameters - is very significant 
(at more than 99.9%). The temperature in this case decreases to 11.2 ±0.8 keV. The photon 
index is 1.6 ± 0.3 and the non-thermal fiux in the 20-100 keV energy band is (0.5 ± 0.2) Fq. 
In the context of statistical significance of raw counting rates, at energies above 20 keV - 
where the power law model dominates - the excess (over the value predicted by the thermal 
model) corresponds to ~ 0.8a and 2.8ct for HEXTE Cluster A and B, respectively. For 
the purpose of illustration, we show the unfolded spectrum of the data fitted to a single- 
temperature thermal model in Figure 2, while Figure 3 shows the confidence levels of the two 
fit parameters, photon index vs. the fiux of the non-thermal component, for our preferred 
two-component model. 

Alternatively, this secondary, hot component can be modeled via another thermal 
plasma spectral component. Using this parametrization (and assuming that only the tem- 
perature and normalization differ from the lower T component), we inferred its temperature 
to be 50 (> 30) keV. The 20 - 100 keV fiux of this component is also 0.5 ± 0.2 x Fq. Now the 
lower T component has a temperature of 10.1 ±0.9 keV. While we cannot clearly exclude this 
interpretation on the quality-of-fit grounds (x^ is 1471 for 1506 d.o.f.), we argue in Sec. 4 
that the power law spectral shape provides for a more viable interpretation of the secondary 
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Fig. 2. — Unfolded spectrum of the cluster 1E0657-56, where the XMM-Newton and RXTE 
data were fitted with a single-temperature thermal (mekal) model, plotted here as the solid- 
line histogram. For the purpose of clarity of the plot, prior to plotting, the XMM-Newton 
data were scaled by a cross-normalization factor, close to unity, to agree with the RXTE 
data. 
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Fig. 3. — Confidence regions (corresponding to Xmin + 2.3,4.6, and 9.2) on the best fit 
parameters for the spectral fits to the joint XMM-Newton , RXTE PCA, and HEXTE data 
for cluster 1E0657-56 for the photon index of the non-thermal X-ray component and its 20 
- 100 keV flux (in units of 10~^^ erg cm~^ s~^). See the text for details. 



4. SUMMARY AND DISCUSSION 

Radio observations of diffuse emission from ICM show presence of non-thermal activity 
in many clusters, especially in those with high SXR emitting gas temperatures (and lumi- 
nosity) and showing recent merger activity. A list of such clusters is shown in Table 1. 
Detections of HXR emission exceeding the levels expected from the thermal gas have been 
reported by two groups using different instruments. This enables a quantitative investigation 
of the nature and origin of the non-thermal activity. The radio observations alone indicate 
presence of extreme relativistic electrons (7 ~ 10'*). Assuming cquipartition between the 
relativistic electrons and the magnetic field we estimate a volume-averaged magnetic field 
value in the range of 0.5 to 2 jiG (see Table 1), in rough agreement with (hne-of-sight av- 
eraged) field strengths deduced from Faraday rotation measurements. In the case of Coma 
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the magnetic field deduced from Faraday rotation of 2 — 3/xG, implies conditions far from 
equipartition if we assume a homogeneous source; the field energy density is about 25 times 
larger than that of the electrons. This discrepancy can be resolved by an inhomogeneous 
model; e.g. the two phase model proposed by Brunetti et al. (2001) that implies a magnetic 
field profile, or if there was 25 times more energy in nonthermal protons than electrons. 

If the HXR radiation is due to IC scattering of CMB photons by the the same electrons 
that produce the synchrotron radio emission, which we have argued to be the simplest 
scenario, such high magnetic files would imply IC fiuxes in the 20 — 100 keV range of about 
10~^Fo while both Coma and A2256 show fiuxes higher than Fq. Production of HXRs at this 
level require magnetic fields of about 0.1 /xG (or partition parameter ( ~ 10'^) which raises 
the required energy of the non-thermal electrons {S^ C^^'^). For example for Coma and 
A2256 the observed HXR fluxes require £e — 2x 10~^^ and 1.5 x 10~^^ erg cm"^, respectively, 
which are comparable to the energy density of the CMB (4.6 x 10"^^ erg cm~^). Assuming 
similar conditions for other clusters, several years ago we identified 1E0657-56 and A2219 
as possible candidate for having detectable HXR fluxes. Here, we have described these 
observations of 1E0657-56 and the results of our spectral fitting to these data alone and to 
the RXTE data combined with the SXR data from XMM-Newton. In the future, attempts 
should be made to observe other clusters in Table 1, and any newly discovered cluster with a 
relatively high temperature and redshift, an ample sign of substructure and/or mergers, and, 
of course, with a strong diffuse radio emission. A low value of magnetic field indicated by 
Faraday rotation measurements or estimated via equipartition argument will also be good 
justification for HXR observations. 

We have shown that a thermal component plus a non-thermal tail described as a power 
law with spectral index F = 1.6 containing ~ 0.5Fo of the 20 — 100 keV flux provides an 
acceptable flt. This flux is fortuitously close to the value of 0.52 Fq we estimate by scahng the 
Coma and A2256 observations. (Note that the equipartition {i.e.( — 1) magnetic fleld and 
flux for this cluster would be ~ 1 /xG and 0.02 Fq.) This agreement is somewhat encouraging 
so that we believe that the chances are even better for detection of HXR emission from Abell 
2219, a good candidate at a moderate redshift, and a recently discovered cluster, MACJ7017, 
with even higher z than 1E0657-56. As is the case for Coma and Abell 2256, these clusters 
also require high nonthermal electron energy densities; £e = [1.0, 1.8,8.5] x 10~^^ erg cm~^ 
for Abell 2219, 1E0657-56 and MACSJ0717, respectively These are even closer to the CMB 
densities at their respectively higher redshifts. In general the required electron energies are 
about 1 to 10 % of the energy density of the thermal gas of few times 10~^^ erg cm~^ which 
is also comparable to the gravitational potential energy density of the cluster. This means 
a very efficient acceleration process. 
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Regarding the spectral fits to 1E0657-56, as shown in Table 2, we also obtain an accept- 
able fit with a double temperature thermal model with a second component with kT — 50 
keV and a normalization implying a volume emission measure {EM) for this component 
which is a significant fraction (~ 10%) of that of the lower temperature component. Indeed 
Markevitch (2005) finds evidence for such a high T component behind the shock near the 
bullet based on the recent Chandra Observatory data. However, this has a much smaller EM 
than is required for fitting the data with a two temperature model. In fact a high T compo- 
nent with such a large EM would have been easily detected by Markevitch (2005). Clearly 
more and deeper observations of the clusters in Table 1 can clarify this situation considerably. 
In particular, deep observations with imaging HXR instruments - such as those with the 
proposed NuSTAR satellite - which can provide information about the spatial distribution 
of the HXRs would be very valuable. 

In spite of the uncertainties about the exact character of the observed radiation spectra 
and the emission mechanism, and in spite of the meagerness of the data, the acceleration 
mechanism of electrons can be constrained significantly. The lifetimes of electrons with 
energies in the range 200 keV < E < 200 GeV are longer than the crossing time, Tcross ~ 
3 X 10^ yr. Therefore, these electrons will escape the cluster and radiate most of their energy 
outside of the cluster unless there exists some scattering agent with a mean free path Agcat ~ 1 
kpc to trap the electrons in the ICM for at least a timescale of Tesc = {R/ Ascat)^cross ~ 3 x 10^ 
yr, for cluster size R ~ IMpc. Turbulence can be this agent and can play a role in stochastic 
acceleration directly, or indirectly in acceleration by shocks, presumably arising from merger 
events. Several lines of argument point to an ICM which is highly turbulent, and there has 
been considerable discussion of these aspects of the problem in the recent literature (see, 
e.g., Cassano & Brunetti 2005). These aspect are explored in POl, but the upshot of which 
is that we require injection of high energy electrons, presumably from galaxies and AGNs, 
and that the acceleration process is episodic on time scale of ~ 10^ to 10^ yr. Whether 
the electrons are injected directly by past and current AGN activity, or are produced by the 
interaction with the thermal gas of cosmic ray protons escaping the galaxies would be difficult 
to determine now. There are however constrains on both these scenarios (see e.g. Blasi 2003). 

The work was supported by NASA's ADP grant NNG04GA66G-NCX, by Chandra 
grants G01-2113X and G04-5125X from NASA via Smithsonian Astrophysical Observa- 
tory, and by the Department of Energy contract to SLAC no. DE-AC3-76SF00515. We 
acknowledge the help in reducing the XMM-Newton data by Mr. Karl Andersson, and Dr. 
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